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Disordered two-dimensional superconductors: 
roles of temperature and interaction strength 

Felipe Mondaini/ Thereza Paiva,^ Raimundo R. dos Santos,^ and R. T. Scalettar,^ 

^Instituto de Fisica, Universidade Federal do Rio de Janeiro, 

Caixa Postal 68528, 21941-972 Rio de Janeiro RJ, Brazil 
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(Dated: July 24, 2008) 

We have considered the half-filled disordered attractive Hubbard model on a square lattice, in 
which the on-site attraction is switched off on a fraction / of sites, while keeping a finite U on the 
remaining ones. Through Quantum Monte Carlo (QMC) simulations for several values of / and U , 
and for system sizes ranging from 8 x 8 to 16 x 16, we have calculated the configurational averages 
of the equal-time pair structure factor Ps, and, for a more restricted set of variables, the helicity 
modulus, ps, as functions of temperature. Two finite-size scaling ansatze for Ps have been used, one 
for zero-temperature and the other for finite temperatures. We have found that the system sustains 
superconductivity in the ground state up to a critical impurity concentration, fc, which increases 
with [/, at least up to [/ = 4 (in units of the hopping energy). Also, the normalized zero-temperature 
gap as a function of / shows a maximum near / ~ 0.07, for 2 < [/ < 6. Analyses of the helicity 
modulus and of the pair structure factor led to the determination of the critical temperature as 
a function of /, for [/ = 3, 4 and 6: they also show maxima near / ~ 0.07, with the highest Tc 
increasing with U in this range. We argue that, overall, the observed behavior results from both 
the breakdown of CDW-superconductivity degeneracy and the fact that free sites tend to "push" 
electrons towards attractive sites, the latter effect being more drastic at weak couplings. 

PACS numbers: 74.20.-z 74.81.-g 74.25. Dw 74.78.-w 
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I. INTRODUCTION 



The interplay between impurities and superconducting 
pairing has been a challenging problem for some time^^ It 
has been argued^ that as long as electronic states remain 
extended in the presence of weak disorder, superconduc- 
tivity should not be affected; as disorder increases, how- 
ever, superconductivity must eventually be suppressed. 
The consequences of this competition are especially in- 
teresting in two dimensions, since the superconducting 
transition belongs to the a;y-model universality class (the 
Kosterlitz-Thouless transition to a state without long 
range order), while non- interacting electrons become lo- 
calized in the presence of any amount of disorder. In- 
deed, by varying the thickness of thin films of Bi and 
Pb deposited on Ge substrates;^ the sheet resistance R^ 
shows insulating behavior (i.e., it increases with decreas- 
ing temperature T) for the thinner samples, and super- 
conducting behavior for the thicker samples; the separa- 
trix between these two regimes extrapolates to a quantum 
critical point as T — > 0,^ and one of the unresolved issues 
is whether or not the threshold of i?n , R* , is universal, 
i.e., R* ^RQ = h/ie^. 

From the theoretical point of view, the two- 
dimensional behavior has been examined in a variety 
of ways. Some have exploited a bosonic description of 
Cooper pairs, ^J^ according to which electrons near the 
Fermi surface are paired and localization is driven by 
Coulomb repulsion amongst pairs. On the other hand, 
one may envisage a fermionic mechanism: disorder en- 
hances Coulomb repulsion amongst electrons, thus de- 
creasing the effective screening (due to electron-phonon 



interactions, in conventional superconductors) which in 
turn leads to the disappearance of Cooper pairs. ^ Inter- 
mediate pictures have also been proposed^, within a phe- 
nomenological theory to address the issue of universality 
of i?*. 

Since none of these approaches have succeeded in fully 
explaining experimental data, alternative routes should 
be sought. One possibility is to study simplified micro- 
scopic fermionic models in which disorder is incorporated 
in a fundamental and unbiased way. However, not much 
is known about models in which disorder is present in the 
pairing interaction. With the purpose of bridging this 
gap, and due to the fact that disorder is more readily 
dealt with in real-space, here we consider the disordered 
attractive Hubbard model, whose Hamiltonian reads. 



{rr')cr rcr 



(1) 



where cl^{c^^) are fermion creation (destruction) oper- 
ators at site r with spin a, and rira = cl^c^^, and H.c. 
stands for Hermitian conjugate of the previous term. The 
kinetic energy lattice sum (rr') is over nearest- neighbor 
sites on a two-dimensional square lattice, and fi is the 
chemical potential; the hopping integral sets the energy 
scale, so we take i = 1 throughout this paper. The on- 
site attraction U{r) is chosen to take on the two values 
C/(r) — and U{r) — U with probabilities / and 1 — / 
respectively; note that C/ > corresponds to attraction, 
according to our definition of the on-site term in ([T]) . A 
review of the homogeneous model can be found in Ref. |8|, 



while Ref. H deals with a recent extension of the model 
to describe non-random layered superconductors, as the 
borocarbides. 

The above model mimics the thin films of Bi and Pb 
referred to above, in the sense that the inverse film thick- 
ness tracks the concentration, 1 — /, of attractive sitesJ^ 
Additionally, it also describes the effects of negative- 
U centers, which are thought to be relevant to high- 
temperature superconductivity in the cupratesi^iiS This 
model has been studied at mean-field level ^^li^iii^ and the 
main results, for a given electronic density, can be sum- 
marized as follows: (i) superconductivity in the ground 
state is destroyed for impurity concentrations above fc] 
(ii) /c decreases as U increases [very simple heuristic 
argumentsii lead to /c = 1 — {U/WY in two dimen- 
sions, with W = 9>t being the bandwidth]; and (iii) 
Tc{f)/Tc{Q) is a concave function of /, which vanishes at 
fc{U). Though mean- field approximations are useful as 
a first approach to the problem, one should be extremely 
cautious about their predictions for two dimensional sys- 
tems. For instance, for the pure system at half filling, 
the degeneracy of charge-density wave and superconduct- 
ing order leads to an effective three-component order 
parameter, thus suppressing the critical temperature to 
by virtue of the Mermin- Wagner theoremjii 



zeroiii^ii^ 



mean-field approaches are unable to detect this feature, 
and should therefore lead to unreliable results close to 
half filling. Indeed, recent Quantum Monte Carlo (QMC) 
simulations have predicted that a small amount of disor- 
der at half filling initially enhances superconductivity;^^ 
this was attributed to the impurity-induced breakdown 
of the above-mentioned degeneracy. In view of all this, a 
more thorough investigation of the model at half filling 
is clearly in order. Here we report on Quantum Monte 
Carlo studies of the dependence of fc with U as well as 
of the dependence of T^ with /, for different values of f/; 
we recall that Ref. 'Tsl was restricted to f/ = 4 and T — Q 
only. As we will see, our predictions are very different 
from those of mean-field approaches. 

The paper is organized as follows. In Sec. [Til we out- 
line the QMC method, and discuss the quantities used 
to locate the superconducting transitions. In Sec. IIIII 
we present a finite-size scaling (FSS) analysis of data for 
superconducting correlations in the ground state, from 
which we extract the behavior of fc with U . In Sec. IIVI 
we perform finite-temperature FSS analyses of data for 
superconducting correlations, as well as analyses of the 
superfiuid density, to obtain Tc{f) for different values of 
U. And, finally, section FVl summarizes our findings. 



II. THE COMPUTATIONAL APPROACH 

We use the determinant QMC methodi^i^aiSii^i^^ to 
investigate the ground-state as well as finite-temperature 
properties of the model. In this approach, the imaginary- 
time interval (0; /3) is discretized into M slices separated 
by the interval At and a path integral expression is writ- 



ten down for the partition function Z. The electron- 
electron interactions are decoupled by the introduction of 
a Hubbard-Stratonovich fieldi^ The fermion degrees of 
freedom can then be integrated out analytically, leaving 
an expression for Z which involves an integral over the 
Hubbard-Stratonovich field, with an integrand which is 
the product of two determinants of matrices of dimension 
the system size. We perform the integral stochastically. 
In the case of the attractive Hubbard model considered 
here, the traces over the spin up and spin down electrons 
are given by the determinant of the same matrix, the in- 
tegrand is a perfect square, and hence there is no sign 
problem»2LS 

In order to study the physics at a particular lattice size 
L X L and value of /, we randomly choose fL^ sites and 
set J7 = on those sites. We typically use 30-50 such 
realizations to average over the different disorder config- 
urations. If fL^ is not an integer, we average over the 
two adjacent integer values, with appropriate weights. 
For each disorder configuration, observables are evalu- 
ated as the appropriate combinations of Green's func- 
tions, which are given by matrix elements of the inverse 
of the matrix appearing as the integrand ii2i2iiS^i2^ The 
average over different disorder configurations then yields 
the quantities of interest. Systematic errors in the calcu- 
lated quantities, associated with our choice of At for the 
discretization of /3, are typically smaller than both the er- 
ror bars associated with the statistical fiuctuations for a 
single disorder realization, and the error bars associated 
with sample-to-sample variations. 

As discussed previously^^ a useful quantity to locate 
the transition is the configurationally-averaged equal- 
time pairing structure factor. 



P, 



Er(r) 



(2) 



where [. . .] denotes average over disorder configurations 
(thus restoring translational invariance) , and the pairing 
correlation function is 

r(r) = (A(i)At(i + r)+H.c.) , (3) 

where (...) denotes ensemble average, with 

A(r)=CrXCrT. (4) 

The different scaling behaviors of Pg, in the ground 
state and at finite temperatures, will be discussed in sub- 
sequent sections. 

Further, current-current correlations probe the super- 
fiuid weight and provide an alternative way to detect the 
destruction of superconductivity. ^4, ^g define 



A,,(r,T) = 0,(r,T)j,(0,0)), 



(5) 



where 



jxir,T) =e 



Ht 



iiV(c, 



/ ^ V^r+^.fj r.fT r,ij r-\-x,(T ) 



^-Hr 



(6) 




FIG. 1: (Color online) Configurationally-averaged equal-time 
pair structure factor Ps as a function of the inverse tempera- 
ture (3, for different lattice sizes at half filling and concentra- 
tion of free sites / = 1/16. 



and the Fourier transform in space and imaginary time, 

r/3 



1 /■'^ 

J^xxiq, ^«) = ITT X! / '^'^^ 



iq-r —iuJnT 



^.Ar,r), (7) 



where Ns is the number of lattice sites, and a;„ — 2mr/p. 
The longitudinal part of the current-current correla- 
tion function satisfies the f-sum rule, which relates its 
value to the kinetic energy K^, 



A =\imq^^o^xxiqx,qy =0,^n = 0) 



A^ 



K^ 



(8) 



(9) 



^r,cr^r+ 



-J). Mean- 



where K^ = {-tEai4+x,a'^r,o 

while, in the superconducting state the transverse part, 

A'^ = lim^^^o ^xx{qx = 0, qy,ujn = 0) , (10) 

can differ from the longitudinal part, the difference being 
the superfluid stiffness Dg, 



Ds/tt - [A^^ - A^] = [K^ - A" 



(11) 



Thus the current-current correlations provide an alterna- 
tive, complementary method to the equal time pair cor- 
relations for looking at the superconducting transition. 



III. GROUND STATE PROPERTIES: fdU) 

As remarked above, there is no 'sign problem' for the 
attractive Hubbard model, so we can do computations 
at very low temperatures (large /3), as shown in Fig. [T] 
for the data of unsealed Ps- We recall that throughout 
this paper we only consider the case of a half-filled band. 



The error bars result from the dispersion in the average 
values taken over disorder configurations. 

The finite-size scaling behavior of Pg allows us to ex- 
tract quantitative information about the superconduct- 
ing transition. As shown by Huse,^^ the spin-wave cor- 
rection to the pair structure factor in the ground state is 
expected to be inversely proportional to the linear lattice 
size. 



L2 



= |Ao| 



a 
Z' 



(12) 



where Aq is the superconducting gap function at zero 
temperature, and a = a{U, /) is independent of L. 
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FIG. 2: (Color online) Zero-temperature scaling of the 
configurationally-averaged equal-time pair structure factor, 
for U — 2, and for different disorder concentrations, /. For 
each /, the intersection with the vertical axis yields the 
(squared) zero-temperature gap; see Eq. (|12|l . 

In Figs. [2H5] we plot the T ^ extrapolated values 
of Ps/ L^ versus 1/ L, for linear lattice sizes ranging fronr 
L = 8 to L = 16, and C/ = 2, 2.5, 3 and 6 (data for 
U — A can be found in Fig. 4 of Ref. il^. According to 
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FIG. 3: (Color online) Same as Fig. El but for U = 2.5. 
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FIG. 4: (Color online) Same as Fig. [H but for U = 3. 
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FIG. 5: (Color online) Same as Fig. but for [/ = i 



Eq. P^ . each intercept with the vertical axis provides 
an estimate for Aq (the square of the zero-temperature 
gap) for the values of U and / considered. For the pure 
system (/ = 0), Aq is pfotted as a function of U in Fig. 
[6l the observed increase of Aq (at least up to C/ ~ 6), is 
due to an increase in the average site double occupancy. 
Indeed, Fig. [7] shows the double occupancy on attractive 
sites, dA = ("-T^-i ^ 1/4) > as a function of U for both the 
pure system and for several disordered configurations. 
The overall behavior is an increase in cIa with [/; and, 
for a fixed U, this double occupancy increases with /, as 
it can be seen from its strong coupling limit, 



K"i-1/4) = J^, 



for U ^ oo. 



(13) 



In order to compare the effects of disorder for differ- 
ent attraction intensities, for each U we normalize the 
zero-temperature gap. A, by their respective pure sys- 
tem values, Aq; the result is displayed in Fig. [51 For 



FIG. 6: (Color online) Zero-temperature gap for the pure 
system, obtained from the f — data, together with Eq. 
(|12p . The full line is a guide to the eye. 
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FIG. 7: (Color online) Average double occupancy on attrac- 
tive sites as a function of [/, for the pure case and for several 
disordered configurations. Data are for 14 x 14 lattices and for 
/3 = 18. Data for U — 100 correspond to the strong couphng 
limit, Eq.[T3] 



2.5 < [/ < 6, the normalized gaps initially increase with 
disorder, reaching maxima [A/AoJ^jj^^^ around / ~ 0.07. 
It should be noticed that [A/Aq]^^^, in turn, does not 
behave monotonically as a function of U, but displays 
a maximum for U — 3, amongst the values of U exam- 
ined. Another crucial information extracted from Fig. [5] 
is that a maximum of A/Aq is absent for [/ = 6. This dif- 
ferent behavior for larger U therefore indicates that the 
disorder-induced breakdown of CDW-supcrconductivity 
degeneracy is not the only mechanism at play: for smaller 
U, the presence of free sites contributes to a decrease in 
the single occupancy, by 'pushing' electrons to the attrac- 
tive sites. For larger U, the pairs are so tightly bound 
that the relative weight of single occupancy is smaller, 
and disorder has hardly any effect on forcing the electrons 
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FIG. 8: (Color online) Normalized zero-temperature gaps as 
functions of impurity concentration, for different values of U. 
Full lines are guides to the eye. 
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FIG. 9: (Color online) Critical impurity concentration as a 
function of U, obtained from Fig. [S] The full line is a guide 
to the eye. 



to occupy the attractive sites. Figure [7] indeed shows that 
the percentual enhancement in double occupancy due to 
disorder is larger ior U ^ A than for f7 ~ 6. As disorder 
increases, the presence of free sites strongly disturbs pair 
coherence, and the gap decreases. The initial increase 
with disorder has also been predicted for anisotropic su- 
perconductors with mcsoscopic phase separation 1^^*^ 

We can obtain the dependence of fc with U by ex- 
trapolating the data for A/Ag to zero. The intercept 
with the horizontal axis of each of the curves in Fig. [5] 
provides fc for the corresponding U, and the result is dis- 
played in Fig. m the error bars reflect the uncertainties 
in the extrapolations of A/Aq to zero in Fig. [51 It should 
be noted that fc initially (i.e., for U < 2) increases very 



slowly with U, which should be attributed to the fact that 
the pairs are not so strongly bound for small C/, so that a 
small amount of free sites destroys phase coherence. As 
U increases, the pairs become more tightly bound, and a 
larger amount of disorder can be sustained before the free 
sites switch their roles, from pushing electrons onto at- 
tractive sites to that of destroying phase coherence. This 
behavior is in disagreement with mean-field predictions, 
according to which fc should decrease with Uf^ How- 
ever, this switching of roles played by the free sites is an 
effect too subtle to be picked up by approaches, such as 
mean-field ones, which do not incorporate fluctuations in 
a fundamental way. For U > A, fc appears to be de- 
creasing with U, and the agreement with the mean-fleld 
approach would set in. 

We recall that the strong-coupling pure attractive 
Hubbard model (at half filling) can be mapped onto 
an isotropic Heisenberg mode L^^i^^ Nonetheless, the fact 
that in the disordered case, fc displays such strong depen- 
dence with [/ is a clear indication that the mechanisms 
of superconductivity suppression by impurities are very 
different from those occurring in diluted magnetic insu- 
lators, which are driven by classical percolation.>22 That 
is, if one is interested in singling out the geometrical as- 
pects of impure superconductors, a model of correlated 
dilution should be more appropriate. 



IV. FINITE-TEMPERATURE PROPERTIES: 

Tc{f,U) 

Another consequence of the two-component nature of 
the order parameter is that at finite temperatures the 
superconducting-normal phase transition for the pure 
system belongs to the Kosterlitz-Thouless universality 
class. And as such, for < T < Tc, one expects that, 
asymptotically. 



r(r) 



"'J(T) 



(14) 



where T{r) is defined by Eq. ([3|), and ri{T) increases 
monotonically between 7^(0) = and 'ri{Tc) — l/iM^ 
The finite-size scaling behavior of Ps is therefore obtained 
upon integration of T{r) over a two-dimensional system 
of linear dimension L. One then has^^ 



P, = L^-'i^T^^FiL/O, i > 1, T^ r+. 



(15) 



where F{z) is a finite-size scaling function of the variable 
z = L/^, with 



^~exp 



A 



[T - T,)i/2 



(16) 



where A is a constant; in the thermod yna mic limit, one 
recovers Ps ^ £j^^. As discussed in Ref.|la, we can obtain 
estimates of Tc by plotting L~''^^Ps{L, (3) as functions of 
/3, for different L, and by looking for intersections/merges 
of curves for consecutive values of L. This procedure was 



supported by independent estimates of the critical tem- 
perature through calculations of the superfluid stiffness, 
Ds, and using the universal jump at Tc, see Ref. [la for 
details. 

The general aspects of the universality class of the 
superconducting transition should remain valid in the 
presence of disorder, since one still deals with a two- 
component order parameter. Further, numerical evi- 
dence has been gathered for the xy-vaoAei^ showing 
that "qlTc) — 1/4 even in the presence of disorder; this 
is in agreement with the Harris criterion, which essen- 
tially states that disorder is irrelevant (in the renormal- 
ization group sense) if the specific heat exponent, a, is 
positive.^ In view of this, our data analyses for the finite- 
temperature transitions can follow along the same lines 
as those for the pure system, ^^ with both Ps and Da 
now being understood as the configurationally-averaged 
equal-time pair correlation function and superfluid stiff- 
ness, respectively. 

Let us first consider the helicity modulus (HM)f2i 
which is given by 



Ps = 



Ds 

47re2' 



(17) 



where Dg is defined in Eq. (jlip , and we take e = 1 in our 
units. At the KT transition, the following universal-jump 
relation involving the helicity modulus holds i^ 






(18) 



where p^ is the value of the helicity modulus just below 
the critical temperature. Thus, on universality grounds 
we may assume the same holds for a configurationally 
averaged ps, and we can obtain Tc by plotting Ps{T), 
and looking for the intercept with 2T/tt. Figure [TUl shows 
data for U = 4, and we see that the intercepts occur at 
Tc = 0.125±0.015, 0.12±0.01, 0.105±0.008, and O.OSOi 
0.015, for / = 1/16, 2/16, 3/16 and 4/16, respectively. 
Similarly to the pure case,^^ we have found here that 
finite-size effects are not too drastic, leading to essentially 
the same estimates for T^. These estimates appear as 
empty circles in Fig. 1111 

In Fig. [T2] we show the scaled configurationally- 
averaged equal-time pair structure factor as a function 
of the inverse temperature /?, for U = A and different 
system sizes, for a given concentration of disorder. For 
usual second-order phase transitions, similar curves for 
two successive linear lattice sizes should cross at a sin- 
gle point, thus leading to estimates for critical inverse 
temperatures. For Kosterlitz-Thouless (KT) transitions, 
on the other hand, curves for different (but sufficiently 
large) lattice sizes should merge above a certain f3c^ In 
the present case, we estimate Pc as the smallest value 
for which the curves for the smallest size superimpose, 
within error bars, with the one for the largest size; this 
ensures that the error bars for data corresponding to in- 
termediate sizes will also superimpose. Thus, applying 
this criterion to the data in Fig. [T^] yields Pc = 8± 1; this 
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FIG. 10: (Color online) Configurationally-averaged helicity 
modulus as a function of temperature, for f/ = 4 and different 
impurity concentrations /, for a 12 x 12 lattice. In each panel, 
the straight line corresponds to 2T/n. 
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FIG. 11: (Color online) Critical temperature (in units of the 
bandwidth) for superconductivity, as a function of impurity 
concentration, /, for different values of the on-site attraction 
U. Full symbols have been obtained through the scaling of the 
pairing structure factor (Ps), whereas the empty circles cor- 
respond to data obtained through the helicity modulus (ps). 
Full lines are guides to the eye. 
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FIG. 12: (Color online) Scaled configurationally-averaged 
equal-time pair structure factor as a function of the inverse 
temperature /3, for different lattice sizes {L x L), at half fill- 
ing, concentration of free sites / = 2/16, and U — 4. The 
curves are labelled by the linear lattice size L, and full lines 
are guides to the eye. 



procedure is systematically repeated for other values of /, 
and we obtain the Tc{f) data for [/ = 4 shown in Fig.fTT] 
The estimates thus obtained are in excellent agreement 
with those obtained from the HM, thus adding credence 
to our merging criterion. 

Given the fact that the calculations of configurationally 
averaged helicity moduli are very consuming in terms 
of computer time (for a given disorder configuration 
the CPU time is increased significantly due to the r- 
integration, and one performs averages over typically 50 
disorder configurations), for other values of U we only 
use data for Pg to estimate Tc{U, /). From Figs. [T51and 
[ni (which yield /3c = 12 ± 1 and (3c = 6 ± 1, respectively), 
as well as from similar ones for other values of /, we ob- 
tain the estimates for Tc{f) for [/ = 3 and U — 6 shown 
in Fig. [m 

Several interesting physical features emerge from Fig. 
[m Firstly, Tc{0) = as a result of the CDW- 
superconductivity degeneracy at half filling; however, any 
finite amount of disorder breaks this degeneracy and Tc 
rises. Secondly, in all curves, Tc displays a maximum at 
some /max, as a result of the interplay between the above- 
mentioned degeneracy and the behavior of the smallest 
energy scale, A. Thirdly, we expect Tc — ^ at fc{U), 
since above /c, superconductivity cannot be sustained 
even in the ground state. And, finally, near /c, Tc dis- 
plays the convex shape observed in experiments^^ the 
steepness of the decrease in Tc{f) can therefore be used 
to fit an effective U by experimental data. As a final 
comment, one should have in mind that some of these 
results should change drastically as the system is doped 
away from half filling. Indeed, since in this case CDW- 
superconductivity degeneracy is already broken in the 
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FIG. 13: (Color online) Same as Fig. [H but for f/ = 3. 
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FIG. 14: (Color online) Same as Fig. [TH but for f/ = 6. 



pure system, Tc should display a monotonic decrease with 
/ for a given [/; nonetheless, we can still expect Tc to be 
a convex function of /. 



V. CONCLUSIONS 

We have addressed the issue of disorder in two- 
dimensional superconductors. To this end, we have con- 
sidered a simple model, namely the attractive Hubbard 
model, in which the on-site attraction is switched off on a 
fraction / of sites, while keeping a finite U on the remain- 
ing ones; the model is defined in such a way that [/ > 
in the attractive case [see Eq. ^\. Through Quantum 
Monte Carlo simulations for typically 50 disorder config- 
urations, we have calculated the configurational averages 
of the equal-time pair structure factor and, for [/ = 4, 
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the helicity modulus, as functions of temperature; there 
is no minus-sign problem in the attractive case. The 
continuous 0{2) symmetry of the superconducting or- 
der parameter allows us to use a spin-wave-like finite- 
size scaling form for the ground state behavior, from 
which the zero-temperature gap was calculated; at finite 
temperatures, the usual finite-size scaling form for the 
Kosterlitz-Thouless transition was used to calculate the 
critical temperature, which was checked for consistency 
against data for the helicity modulus. Our numerical 
data are consistent with the following findings: (i) Su- 
perconductivity in the ground state is destroyed above 
an impurity concentration, fc, (ii) At least up to C/ = 4, 
this critical concentration increases with increasing U, 
slowly for U < 2.5, and then fast up to U ^ 4; this be- 
havior does not agree with mean-field predictions, due to 
the important role played by fiuctuations, not included 
in the latter approach. The error bars prevent us from 
ascertaining that fc decreases with U above U — A, but it 
may be that the mean field behavior is recovered in this 
regime, since fluctuations should become less important 
for large U. At any rate, the transition at zero tempera- 
ture is not driven by purely geometrical aspects, such as 
in dilute insulating magnets, (iii) In the range between 
U = 2.5 and U <6, the normalized zero-temperature gap 
initially (i.e., small disorder) increases with disorder, as a 
result of both the breakdown of CDW-superconductivity 



degeneracy and the fact that free sites "push" the elec- 
trons towards attractive sites; and, (iv) near the critical 
concentration of defects beyond its maximum value, Tc 
is a convex function of /, as observed in experiments; 

Overall, we conclude that the random attractive Hub- 
bard model is a promising working ground to investigate 
the interplay between impurities and pairing. By tuning 
two variables at half filling, namely the impurity concen- 
tration and the pairing potential, we have found instances 
in which small disorder either hardly affects supercon- 
ductivity or enhances it. It should therefore be of in- 
terest to check whether these features remain valid away 
from half filling. Further, the present model can be used, 
with suitable changes, to investigate other disordered 
BCS superconductors, such as three-dimensional carbon- 
substituted MgB2,^^ and MgB2/MgO superstructures!^ 
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